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Trajectory of a test particle or a photon around a general spherical black hole is studied and 
bending of light trajectory is investigated. Pseudo-Newtonian gravitational potential describing the 
gravitational field of the black hole is determined and is compared with the related effective potential 
for test particle motion. As an example, results are presented for Reissner-Nordstrom black hole. 

I. INTRODUCTION 

Motion of a test particle around a black hole is a very old topic of investigation to know the behavior 
of the gravitational field around black hole(for example [1-3] and any other book on general relativity). 
Such investigation was started long back when within a year of publication of general theory of relativity 
by Einstein(1915)[4-7], Schwarzchild [8] gave a vacuum solution to the Einstein field equations. The 
solution describes the geometry of vacuum space-time outside a spherical massive body and is known 
today as Schwarzchild black hole solution. At present in all text books on general relativity there 
is an exhaustive study of the motion of a test particle around the Schwarzchild black hole (see for 
example [1-3]). In the present work, an investigation of the motion of a test particle is done around a 
general static non-rotating black hole. A general formula for determining bending of light is evaluated 
and is tested for Schwarzchild black hole. Gravitational field outside the black hole is approximated 
by Pseudo-Newtonian(PN) gravitational potential and is compared with the corresponding effective 
potential for test particle motion. Finally all results are verified for Reissner-Nordstrom black hole 
solution. 

The paper is organized as follows: Section II deals with the motion of a massive test particle around 
a general spherically symmetric non-rotating black hole. Also effective potential, energy and condition 
for circular orbit are determined for the black hole. Further some comments are presented point wise. 
In section III trajectory of a photon and bending of light is studied. Pseudo-Newtonian gravitational 
potential is determined and it is compared with effective potential in section IV. As an example all results 
are deduced for Reissner-Nordstrom black hole in section V. The paper ends with a short conclusion in 
section VI. 



II. MOTION OF A TEST PARTICLE AROUND A GENERAL BLACK HOLE 



The line element of a static spherically symmetric space time which describes a black hole can be 
written as, 

d^ = -f(r)dt 2 + ^-+r 2 dni (1) 

where f(r) is at least a c 2 function and it should satisfy the following conditions so that line element 
(eq.l) describes a black hole solution: i) /(r) must have a zero at some positive r (say r^) so that time 
dilation is infinite at r^, ii) The Kretschmann scalar (a = R l i kl Rijki) should be finite at r = but it 
diverges at r = i.e. the space time described by equation (1) has curvature singularity only at r = 0. 
Here 

dflj = d6 2 + sin 2 9d(j) 2 
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is the metric on unit two sphere. Suppose we consider the motion of a test particle of rest mass m around 
the black hole. The corresponding Lagrangian will be 



2 / , \ 2 / m \ 2 



dO 



\dX) ^ JV> \dX) ' ' \dX) " " \dX 
where A is any affine parameter. 

As the Lagrangian has two cyclic co-ordinates t and <f) so the corresponding momenta must be constant. 
This leads to 

E = -™ (3) 
m 



a constant and 



L = — (4) 

m 



is also a constant. 

As the space-time is spherically symmetric so the motion is always confined to a plane which for 
convenience chosen to be the equatorial plane = f . The explicit form of the momentum components 
are 

p° = E 



m 

W) 

- dr 

p =m Tx 

V 9 - (5) 



r z 

Using the above expressions for momentum components in the energy momentum conservation relation 

p"p„ = -m 2 (6) 

we obtain 

(!)'-*-"» (7) 

Where 

V 2 (r) = /(r)(l + ^) (8) 

is called the (square of the) effective potential. 

Now differentiating both sides of equation (7) we have 

d 2 r l dV 2 (r) 
dX 2 ~ ~2 dr 

Also from (5) the momentum in the <j> direction gives 

d(j) L 



(9) 



dX - r 2 (10) 
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So eliminating the affine parameter between (7) and (10) the differential equation of the trajectory of the 
particle in the equatorial plane is given by 



c?r x 2 



dtj>) L 2 



(11) 



Which can be written as 



Where 



dr\ 2 r 4 



<P(r)=E 2 -f(r){l + ^r) (12) 
We can make the following conclusions on the trajectory of the particle: 

• The energy of the particle should not be less than the potential V(r) i.e. for a given E the 
trajectory should be such that the radial range is restricted to those radii for which V is smaller than E. 

• If i/j(r) > for all values of r then the particle comes from infinity and moves directly to the 
origin. This is called the terminating escape orbit. 

• If ijj(r) has one positive zero then the particle either starts from finite distance moves directly to 
the origin (known as terminating bound orbit) or it may move on an escape orbit with a finite impact 
parameter ^. 

• If ip(r) has two positive zeros then we have two possible cases: I. if ip(r) > between the two 
zeros then the trajectory is called periodic bound orbit like planetary orbit or II. if tp(r) < between 
these two zeros then the trajectory is either an escape orbit or a terminating bound orbit. 

• The points where ip{r) = are known as turning points of the trajectory i.e. the value of r which 
satisfies 

e 2 = m (i + (13) 

are turning points and Eq.(13) determines the potential curves. 

• For circular orbit (r = constant) we have from (9) rf ^ r ^ = i.e. circular orbits are possible for 
those radial co ordinates which correspond to maximum (unstable) or minimum(stable) of the potential. 
Thus for circular orbit we must have 

™ = ^ (14) 
f(r) r(r 2 +L 2 ) 1 ' 

III. TRAJECTORY OF A PHOTON: BENDING OF LIGHT 

To determine the photon trajectory we shall proceed as before. Here from the energy momentum 
conservation relation we have 

dr\ 2 „, S L 2 



i.e. Vl{t) 2 = f{r)^r is the (square of the) effective potential. So differentiating both sides we get 

d 2 r l dV 2 (r) 

d\ 2 2 dr { 1 
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Thus the differential path of a light ray is given by 

d<f> 



, =± , (17) 



where b — h. Now for photon circular orbit we have 



dr 

i.e. the radius of the circular orbit satisfies 

rf'(r) = 2/(r) 

One may note that the radius of photon circular orbit is independent of the angular momentum of 
photon. In particular, for a Schwarzchild black hole the radius of photon circular orbit is 3M .For an 
ingoing photon choosing u — £ we have 

where F(u) = /(I). 

Note that if F(u) is a constant then (18) has the solution 

r sm((/) - 4> Q ) = b (19) 

(choosing the constant to be unity) a straight line. This is expected as f(r)=constant means the space 
time is minkowskian (having no gravitational effect) and photon trajectory will be straight line. Further 
we see that at large distance (small u) the gravitational field due to the black hole will be negligible so 
we may expand F(u) in a power series of u i.e. 

F(u) = 1 + cm + c 2 u 2 + ■■■ (20) 

So keeping up to first order in u we have from (18) 

T - i 1 (21) 

Let y = u(l + £ y l ) then u = y (l — ^)+0 (u 2 ). Then the above differential equation becomes 

^ = 4^L+0(«») (22) 

So on integration, 



^ = 0o - j + sin" 1 (by) + Cl ^1 - j/2 (23) 

If we choose 0o as the initial incoming direction of light i.e. <fi — > <f>o as y — > and as in the 
approximation y= \ corresponds to the smallest r that photon can travel then 

*y=i = - J + \ (24) 

So the angle of deflection is 5 — ^ — as the photon comes from infinity to the point of closest approach. 
Hence the total deflection would be S = n — Therefore, considering the straight line path the net 
amount of deflection will be 

A0 = -^i (25) 
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If we take schwarzschild solution then we have c\ = —2M and A<fi = Further, if it so happen that 
ci = then we choose F(u) ~ 1 + c 2 u 2 . Using the transformation 



, ., . c 2 u 2 

y = u[i + —— 



the net amount of deflection will be given by 



A*=^ <») 



IV. PSEUDO-NEWTONIAN GRAVITATIONAL AND EFFECTIVE POTENTIALS 

Based on a general heuristic method([9]) ,the PN gravitational potential can be defined as 

1> = J l ldr (27) 

where, r is the usual radial co-ordinate and l c is the general relativistic specific angular momentum i.e. 
l c (= jjf-) is the ratio of the conserved angular momentum and energy per particle mass, related to the 
circular geodesic in the equatorial plane. 

In newtonian theory, the gravitational potential ip n is given by 



/2 

3 



with l cn , the newtonian angular momentum per mass of the particle moving in the circular orbit. The 
motivation of choosing PN potential (27) is to match the newtonian angular momentum per particle mass 
on a circular orbit with the general relativistic angular momentum. In the present study, the general 
relativistic conserved angular momentum and energy per particle mass for circular orbit are given by 
(from equations (13) and (14)) 

r 3 f'(r) 

JK) (28) 



2/( r ) - rf'(r) 



and 



E C = - ^ (29) 
VWW - rf'(r)\ 



i.e. 



lc ~my~^ (30) 

Hence from (27) the PN gravitational potential is 

* = c -*k ,31) 

where the constant of integration 'c' is determined from the known result of the Schwarzchild black hole 
as follows: (note that c has no physical meaning) 

For Schwarzchild black hole, the well known Paczyhski-Witta gravitational potential([10]) is 

M 

^ = -^~2M (32) 



G 



substituting in eq.(31) we get c = \ and hence the PN gravitational potential for a general spherically 
symmetric black hole described by eq.(l) is 



f(r) 



(33) 



As for static radius (r s ) gravitational potential should be zero so from (33) we have f(r s ) = 1. Hence 
from equations (28), (29) the circular orbit of the test particle exist for radius (r c ) in the range([ll]) 



where r a satisfies 



r a <r c < r s 



rf(r) =' 2/(r) 



i.e. r a is the photon circular orbit. Thus all circular orbits of the test particle are lower bounded by the 
photon circular orbit and are extended upto the static radius. 

Further, from eq.(33) we see that the PN potential diverges at the event horizon (i.e./(r) = 0) reaches 
its maximum value at r = r m (where/' (r m ) = 0) and then decreases for r > r m i.e. the gravitational field 
corresponding to PN potential become repulsive for r > r m . Also if the metric(l) becomes asymptotically 
flat(i.e. f(r) — > 1 as r — > oo) then ijj — > asymptotically ([11], [12]). Moreover, for central gravitational 
fields, if we assume that the motion of the test particle is confined to the equatorial plane then for 
Keplarian motion along the radial direction gives 



1 (dr 



2U ] =E ~ Veff 



(34) 



where e stands for total PN energy per particle mass and v e ff stands for PN effective potential per 
particle mass having explicit form ([13]) 



v eff = i> + 



_L_ 

2r 2 



(35) 



where ip is the PN gravitational potential given by eq. (33) and I is the PN angular momentum per particle 



mass. Thus circular Keplarian orbits are characterized by the extrema of v e ff(i-e. v t? f = 0) and we have 



f = r 3 f(r) 



2[f(r)Y 



which can be written as (using eq.s (28) and (29)) 



The corresponding expression for energy is 



E P . 



r/'(r)-2/(r) 



2(/(r))' 



using E c it can be written as 



1 



1 



(36) 



(37) 



(38) 



(39) 



It is to be noted that angular momentum per particle mass is same for both General Relativity as well 
as for PN effective potential theory. 



We shall now examine the stability of circular orbits studied above. As stability criteria is determined 
by the extrema of the effective potential v e ff, so we have 

dl c n 
IT >0 

or 
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for stable circular orbit and 

9l c n 
IT < 

for unstable circular orbit. Due to identical nature of l c both GR and PN potential theory have same 
criteria for stability. Further, inner and outer marginally stable circular orbits corresponds to cxtrema of 
I 2 and we have 

2r[/'(r)] 2 = /(r)[3/'(r)+r/"(r)] (40) 
Hence for stability one should have 

2r[/'(r)] 2 </(r)[3/'(r)+r/"(r)] (41) 
Finally, from the above expressions (eq.s (36) and (38)) of I 2 and e c we have the following observations: 

• At the event horizon(/(r) = 0) both I 2 and e c diverge while E c and L c exist if [2/(r) — rf'(r)] is 
positive definite. For example in case of Schwarzchild-de Sitter space time([ll]) they become finite while 
for our Reissner- Nordstrom black hole (see next section) they do not exist. 

• At the static radius (/(r) = 1) both I 2 and e c vanish while E c and L c depend on choice of f(r). 

• At the photon circular orbit (rf'(r) = 2/(r)) both L c and E c diverge while e c = \ and l c is finite 
there. 

Thus we conclude that both in PN approach and in rclativistic approach the circular orbits arc 
bounded from above by static radius while bound from below by event horizon in PN approach and that 
by the radius of photon orbit in relativistic approach. This is due to the fact that we do not get photon 
circular orbit in PN approach([ll]). 



V. AN EXAMPLE: REISSNER-NORDSTROM BLACK HOLE 

For Reissner-Nordstrom(R-N) black hole solution we have 

2M q 2 , , 

/M = l-— + ^2 ( 42 ) 

where M being the mass and q the charge of black hole. So for horizon we have f(r) = i.e. 



r± = M± y/M 2 - q 2 (43) 

Black hole solution exists for 

M 2 > q 2 

and we have event horizon at 

rh = r+ 

and 

r c = r_ 

corresponds to black hole Cauchy horizon. When 

M 2 = q 2 

then both the horizons coincide and it is the case of extremal black hole. 
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Figure 1 



Figure 1: The figure shows the variation of relativistic conserved angular momentum L c with respect to 

variation of q and r for the choice M = 1. 



The static radius is given by (f(r) = 1) 



The radius at which f'(r) = Ogives the solution 



2M 



(44) 



One may note that 



M 



r s < r m < r + 



So both the static radius and r m lie inside the event horizon. Hence they have no physical significance 
for R-N black hole. 

Also the radius of photon circular orbits are given by 



1 



3M ± V^M 2 - 8q 2 



(45) 



Relativistic Theory: 

For circular orbit the conserved angular momentum and energy per particle mass are 



T\J Mr — q 2 
yfr 2 - SMr + 2q 2 



(46) 



Variation of L c with respect to the variation of q and r for M — 1 has been shown in figure 1 [variables 
in the figures are in any standard units]. 



E r . = 



2M | q- 
— ~ T" 



y/r 2 - Mir + 2q 2 



(47) 







Figure 2 

Figure 2: Here variation of PN gravitational potential ip with variation of both q and r for M = 1 is shown. 



L c _ y/Mr - q 2 



(48) 



with effective potential 



V, 



eff 



1 



2M 



L 2 



(49) 



Pseudo-Newtonian Theory: 

PN gravitational potential and effective potentials are 



2 [r 2 - 2Mr + q 2 } 



(50) 



The graph of -0 for variation of both q and r with M = 1 is presented in figure 2. Also in figure 3, we 
have shown the dependence of ip by drawing the graphs of ip for three different values of q. 



We// 



2 [r 2 - 2Mr + q 2 } 



2r 2 



(51) 



Also we have 



lr.= 



y/Mr - q 2 
1-2M + 4 



(52) 



The graphical presentation of l c for the variation of both q and r with M = 1 has been shown in figure 4. 
Also graphically a comparative study of L c (given by eq.(46)) and / c (given by eq.(52)) has been presented 
in figure 5(a) and 5(b) for two different values of q. 
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Figure 3 



Figure 3:Here PN gravitational potential ip is drawn for three different values of q. The upper one for q — 0, 
middle one for q = 0.8 and the lower one for q = 1. 




Figure 4 



Figure 4:The figure shows the graphical representation of PN angular momentum per particle mass l c for 

variation of both q and r with M = 1 
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Figure 5(a) 



Figure 5(b) 



Figure 5(a) and 5(b) show the comparative study of the variation of relativistic angular momentum L c (given by 
eq.(46)) and PN angular momentum per particle mass £ c (given by eq.(52)) for q = 0.8 and 0.6 
respectively. Where upper one gives L c and lower one l c 



_ q 4 - Mr 3 - 4Mrq 2 + 4M 2 r 2 
Cc ~ 2[r 2 -2Mr + q 2 } 2 

The marginally stable circular orbits are given by the positive root of the equation 

Mr 3 ~ 6M V + 9Mq 2 r - 4q 4 = (54) 

If the cubic eq has three positive real roots(r m i < r„a < r m 3)then r = r m ^ and r = r m \ are respectively 
the radii of the outer and inner marginally stable circular orbits. 
Finally for stable circular orbit we have 

Mr 3 - 6M 2 r 2 + 9Mq 2 r - Aq 4 > (55) 



VI. CONCLUSION 



In this work we give a general formulation of the trajectory of a test particle (or a photon) around 
any spherically symmetric black hole in four dimensional space time. We also classify the trajectories by 
studying the possible positive zeros of the function 

V(r) = E 2 f(r) (l + 

So once f(r) is given for a given black hole we can immediately tell the trajectories of a test particle 
or a photon around it. In the PN approach physical parameters for circular orbits are evaluated and 
are compared with the corresponding quantities in relativistic treatments. Stability condition for the 
circular orbits are determined and bounds of the marginally stable circular orbits are compared in both 
formalism. As an example we have applied our results for R-N black hole solution. 
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It is to be noted that the above analysis of the trajectory is not restricted to Einstein gravity, it can also 
be applied to any black hole solution in a modified gravity theory. Further our analysis can be extended 
to any higher dimension. For example, the metric ansatz for a n dimensional black hole is written in the 
form 

fir 2 

ds ^-f(r)dt 2 + — +r 2 dnl_ 2 (56) 
f(r) 

Where dCt 2 l _ 2 is the metric on unit n — 2 sphere and is given by 

dnf = d(f> 2 , 



dnf +l = def + sin 2 o t dn 2 ,i > 1 

Then due to spherical symmetry of the space time the motion of a test particle can be restricted to the 
equatorial plane defined by 6i — > 1. As before energy E and the angular momentum L are two 
conserved quantities and the differential equation for the path of the test particle becomes 

where e = or 1 for photon or massive particle. Finally, one may note that throughout our calculations 
we have not used any specific gravity theory. So if we have a black hole solution given by equation (1) 
not only in Einstein gravity but also in any other gravity theory, the above analysis of the trajectory of 
a test particle is valid. Moreover from equations (56) and (57), we may conclude that above analysis of 
particle trajectory can be extended to any dimension of space time. 

For future work, an extension of the above approach to non-spherical systems (particularly axi- 
symmetric) would be interesting. 
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